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Abstract. The semiclassical description of billiard spectra is extended to include the 
diffractive contributions from orbits which are nearly tangent to a concave part of the 
boundary. The leading correction for an unstable isolated orbit is of the same order as 
the standard Gutzwiller expression itself. The importance of the diffraction corrections 
is further emphasized by an estimate which shows that for any large fixed k almost all 
contributing periodic orbits are affected. The theory is tested numerically using the 
annulus and the Sinai billiard. For the Sinai billiard the investigation of the spectral 
density is complemented by an analysis which is based on the scattering approach to 
quantization. The merits of this approach as a tool to investigate refined semiclassical 
theories are discussed and demonstrated. 



PACS numbers: 03.65. Sq, 05.45.+b 
1. Introduction 

Beside the generic orbits which chaotic billiards support, there may exist important 
families of non-generic orbits, which affect the dynamics, and play a prominent role 
when the billiards are quantized semiclassically. The bouncing ball orbits which reflect 
between straight sections of the boundaries are an important example which was studied 
extensively in the past [p], 0] and will not be dealt with here. Another type exists in 
convex billiards (or along convex sections), and it comprises of whispering gallery orbits: 
these are classical trajectories which provide a hierarchy of polygonal approximants to 
the boundary. A detailed study of the role of such orbits in the quantization of convex 
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Figure 1. The quarter Sinai billiard and the attached waveguide. 



and smooth boundaries is given in ||. The particular case of the stadium billiard, to 
which Lazutkin's theory does not apply, was studied in 0, |fl. The whispering gallery 
orbits occupy a narrow strip in phase space, limited on one side by the boundary of 
the phase space domain. In the case of a smooth, convex billiard, every point on the 
phase space boundary is a fixed point of the bounce map. We can therefore consider 
the boundary as a one-parameter family of fixed points, which is the limit of the family 
of whispering gallery orbits. 

In the present paper, we shall focus our attention on concave billiards, where the 
whispering gallery modes do not appear. Instead, concave billiards are characterized 
by the existence of tangent orbits along the concave sections of the boundary. They 
also form a one parameter family, which belongs to the boundary of the phase space 
domain. In this sense they are the counterparts of the whispering gallery orbits in convex 
billiards. In contrast to convex billiards, however, tangency introduces discontinuities in 
the classical map, and therefore the set (of zero measure) of tangent orbits is excluded 
from the classical phase space when the ergodic properties of the billiards are studied. 
Tangency is also responsible for a special kind of bifurcation which can be best illustrated 
by considering the standard Sinai billiard, or rather, its desymmetrized quarter (see 
figure [l]). This will be the example we shall use throughout this article and it can be 
generalized to other chaotic concave billiards in a straightforward manner. The chord 
AB in figure [1] is tangent to the arc of radius Rq and it is a classical trajectory which 
leads from A to B. Bifurcations due to tangency occur when one allows the radius of the 
arc to vary. When the radius is reduced to any value i?_ < Rq, two classical trajectories 
can serve to connect A and B. One of them goes directly, the other reflects specularly 
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from the arc at a point C. As i?_ — > R$ the two trajectories become closer, until they 
coalesce when R_ = Rq. When the radius is increased to values R + > R , both the 
direct and the reflected trajectories from A to B become classically forbidden. That is, 
A and B are mutually shaded from each other by the geometrical shadow cast by the 
arc. 

Tangency affects also the quantum (wave) dynamics in the billiard. Due to the 
finite wave length, the sharp geometrical shadow is replaced by a transition region which 
smoothly interpolates between the strictly illuminated and shaded regions. This is the 
penumbra (Latin: almost shadow) domain. The standard semiclassical quantization 
of billiards, which is restricted to the illuminated domain, expresses the density of 
states d(k) = ^2 n S(k — k n ) (E n = k\ are the eigenvalues) in terms of classical periodic 
orbits via the Gutzwiller trace formula ||. A semiclassical theory for transitions to 
the strictly shaded domain can be formulated in terms of non-classical orbits which are 
allowed to creep along a section of the boundary. This approach was first discussed 
by Keller and co-workers who developed the concept of geometrical diffraction theory 
in a systematic way ||. Vattay et al. |7[ have recently generalized Gutzwiller's trace 
formula by including periodic orbits with creeping sections that give exponentially small 
corrections. In a previous paper || we gave a preliminary account of a semiclassical 
theory which is valid in the penumbra where the diffraction contributions are not 
exponentially small. The purpose of the present paper is twofold - to give a complete 
discussion of the penumbra effects, and to show how they can be scrutinized within the 
scattering approach to quantization [Q, [LT], [□], [T2j . 

Diffraction effects appear as corrections to the leading semiclassical expressions 
which are based on Gutzwiller's trace formula. To identify corrections to the 
contributions from individual periodic orbits, it is advantageous to study the "length 
spectrum" 

/•OO °° 

D(x)= dke ikx d(k) = Ve' 1 " 1 . (1) 

Jo n=l 

Strictly speaking, D(x) is a tempered distribution whose singular support is at lengths 
of periodic manifolds of the classical billiard [T3J. In the Sinai billiard these are the 
continuous families of neutral orbits (bouncing ball manifolds) and isolated unstable 
periodic orbits. In practice, we do not have the complete spectrum at our disposal, and 
we study the truncated length spectrum 

/•oo 00 

D w (x) = dke ikx w(k)d(k) = J2 w ( k n)e iknX ■ (2) 

JO n=l 

Here w(k) is a smooth positive function with a finite support centered at k — k. D w (x) 
is a smooth function with finite spikes where D(x) is singular. The leading semiclassical 
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contributions to (Q) can be evaluated by substituting the Gutzwiller trace formula, 
augmented with the special expressions due to the neutral families. The difference 
between the semiclassical and the exact length spectra gives a measure of the quality 
of the semiclassical approximation. We shall show below that the corrections due to 
diffraction effects are responsible for the largest deviations. The main difficulty of 
this approach is that the length spectrum for a chaotic billiard is rather dense, and 
even with spectral sections containing thousands of eigenvalues, the contributions from 
neighbouring lengths overlap to the extent that a detailed investigation of contributions 
from individual periodic orbits is impossible. However, the scattering approach to 
quantization offers a method to disentangle some of this complexity which is due to 
the proliferation of periodic orbits. 

The scattering approach we shall employ here makes use of an auxiliary scattering 
system which couples the billiard to a wave guide in an appropriate manner [§, [ITj, [TT], [TJ] . 
In a way, this method provides the quantum analogue of the classical Poincare section 
in terms of the scattering matrix, and it has much in common with earlier []14| and later 



15, 16 1 attempts aiming at a similar goal. Since the method is well documented and 



reviewed, we shall only mention its basic ingredients for the particular application to 
the Sinai billiard. 

The scattering problem which we use in the present context is described in figure 
[TJ. We define two systems to which we apply the Krein spectral shift theorem [ 17j . The 
Hamiltonian Ho of the first system is the kinetic energy (—A) in the space of functions 
which satisfy Dirichlet boundary conditions on the parallel channel walls and on the 
section Tg which separates the original Sinai billiard from the wave guide. In the second 
system H the wall at Tg is not present and the functions under consideration have 
to satisfy Dirichlet boundary conditions on the extended billiard boundary displayed 
with a full line in figure (UJ). One can now define the scattering matrix S(E) for any 
energy E = k 2 , and at the energy eigenvalues E n of the billiard the secular equation 
det(J + S(E n )) = is satisfied. Krein's theorem is expressed by the relation 

- hm Im Tr [G(E + ie) - G (E + ie)] = • (3) 

The trace in (H) is taken over the space of continuum eigenstates of H . G and G are the 
Green functions for the two systems (Go is zero inside the billiard domain), and 0(E) 
is the total phase of the scattering matrix S(E), defined by 0(E) = — i logdet S(E). 
Krein's theorem connects the excess density of continuum states due to the introduction 
of the scatterer, with the total phase of the S-matrix. Performing the trace operation 



in the semiclassical approximation (see e.g. [TR |T9||), one finds that the left hand side in 



(0) can be expressed as a sum of a smooth term and contributions from trapped periodic 
orbits - orbits which do not escape in spite of the fact that the billiard is opened. This 
is very important in the context of the semiclassical quantization of billiards by the 
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scattering approach. There, one expresses the spectral density of the closed billiard as 

v ; 2vr dk TT dk ^ n y ' y ' 

In the semiclassical limit, each Tr S n is expressed as a sum of contributions from the 
periodic orbits of the closed billiard which bounce on the section Ts exactly n times. 
However, there may be periodic orbits which do not reflect from Ts at all. They are 
trapped in the open billiard, and their contributions to the level density come from the 
semiclassical expression of j- d< ^^ via @, as explained above |§. d ®j^ provides also 
the leading term in the expression for the smooth part of the spectral density. It contains 



higher order contributions as well, but not necessarily all of them, as is discussed in f20 

Taking the Fourier transform of ([|), we get an expression of the length spectrum 
in terms of the Fourier transforms of and of dTr ^ fc ^ for all positive n. The 

transform of the total phase will provide the length spectrum of orbits which are trapped 
in the open billiard. In our system there exists only one isolated and unstable trapped 
orbit. Another contribution comes from the family of marginally stable orbits which 
bounce perpendicularly between the straight sections of the billiard. This is a very 
sparse set of lengths, and it leaves a lot of space to observe diffractive orbits of various 
kinds. The Fourier transform of dTr ^ fc — provides the length spectrum of orbits which 
bounce n times from the section Tg. In this way we can partition in a systematic way 
the contribution from orbits with different n to the total length spectrum and next to 
leading order effects such as diffraction contributions can be better observed. 

To have an idea of the way how the theory which was described above works in 
practice, let us consider the simple case of the square billiard. Denoting the length of 
the square by a, we get a diagonal scattering matrix 



Si,i>{k) = -exp 




6u> ■ (5) 



The subspace of conducting modes is of dimension A = [— 1 , where the symbol [•] stands 
here for the integer part. In this case 



P-- 2 (6) 



Using the Poisson summation formula and other standard relations, one gets 

J-6(*) = ^-~- + \ + ^f: -M2mak) - i f^-s^mak) . (7) 
Z7T Air 7i 4 2tt ^— ' m lix *— ' m 

m=l m=l 

This is an exact equality which can be interpreted in the following way. The expression 
appearing in the upper line of (0) is the smooth spectral counting function N(k) which 
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consists of three terms only - the area, circumference and corners terms jn|] • The next 
infinite sum is the contribution of the (open) manifold of periodic orbits which are 
parallel to the section IV The last sum is due to to the limiting periodic orbits which 
run along the edge Ts and its counterpart on the other side of the billiard. These two 
limiting orbits are the closure of the manifold mentioned previously. We thus see that 
the oscillatory part of (|7|) consists of contributions which exhaust all possible trapped 
periodic motion in the open billiard. We would like to emphasize again that this result 
is exact. In the sequel, when we treat the more interesting Sinai billiard, we shall obtain 
similar relationships which involve other possible trapped orbits. This will be done, 
however, within the semiclassical approximation and its refinements. 

The paper is organized in the following way. In the next section we shall present 
the semiclassical theory of quantization for billiards which are exterior to a circle taking 
into account the leading diffraction effects both in the penumbra and the deep shadow 
domains. Once this is achieved, we demonstrate the success of the theory using two 
numerical examples. The first example in section ^ will be the annular sector which 
represents an integrable billiard. In the following section ^ we turn to a detailed 
numerical study of the quantized Sinai billiard which is chaotic and therefore necessitates 
the methods of analysis which were outlined above. In section we shall discuss 
the general significance of diffraction corrections for the semiclassical quantization of 
a concave and chaotic billiard. The results of the paper will then be summarized in 
section 0. 



2. Theory of diffraction for dispersing billiards 

In this section, we consider billiards with a domain Q that is exterior to a circle (e.g. the 
Sinai billiard). We find expressions for the contribution to the density of states of 
periodic orbits which are nearly tangent to the circle (either by reflecting with a very 
small angle or by passing very close to the circle). As will be shown, the standard 
semiclassical expressions for contributions of such orbits fail. In addition creeping orbits 
exist due to the concave billiard boundary (the circle in our case). The contribution of 
periodic orbits which have a creeping part was studied by Vattay, Wirzba and Rosenqvist 
0. However, their expressions for creeping orbits also fail, if the orbit is too close to 
tangency (i.e. if the creeping angle is too small). The expressions for periodic orbits 
near tangency, derived in this section, extend the semiclassical description of a general 
billiard exterior to a circle. The methods we will use for nearly tangent periodic orbits 



are adapted from the calculations of Nussenzveig [22] for the problem of scattering off 
a three-dimensional sphere. 

The free Green function satisfies 



(A + k 2 ) G(r, r') = -5(r - r') 



(8) 
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for any r, r' and outgoing boundary conditions at infinity. From Green's theorem one 
obtains that the eigenvalues k n of the billiard are those values for which the boundary 
integral equation 

/dG 
&s' —{r s ,r s >)u{r s >) (9) 

for the normal derivative of the wave function u(r s ) = di/)/dn(r s ) has a solution. The 
billiard boundary is denoted by T, and the normal direction n s in a point on the 
boundary r s is pointing from Q outwards. 

Equation @ is used to obtain a secular equation, from which the density of states 
d(k) may be found by a multiple reflection expansion containing terms as 

x 2N d f , , 9G / \ dG , ,dG / 

Im— — / dsi . . -ds N — r-(ri,r 2 ) • • • ^ {r N - 1 ,r N )——(r N ,r 1 ). (10) 

When the integrals are evaluated in stationary phase approximation, each saddle point 
in one of these terms yields the contribution of a periodic orbit which reflects N times 
off the boundary. This method has often been used as a starting point to derive 
semiclassical theories for billiards, e.g. Alonso and Gaspard employed it for finding higher 
order corrections to the contributions of periodic orbits |23fl. 

For the purpose of considering nearly tangent or creeping orbits at the circle, the 
free Green function is replaced by the circle Green function, which will also be denoted 
by G. It satisfies (§) for r and r' exterior to the circle, and moreover the prescribed 
boundary conditions on the circle. The derivation of (|10D may still be followed and 
yields the contribution of periodic orbits which reflect N times off T, which is now the 
billiard boundary excluding the circle. The reflections off the circle are included in the 
Green function, and so is the diffractive behaviour associated with creeping or nearly 
tangent orbits. The problem is thus reduced to understanding the behaviour of the 
Green function of the circle for different positions of r and r'. 

Although in principle the multiple reflection expansions based on the free and the 
circle Green functions are equivalent and both exact, the latter greatly simplifies the 
mathematical treatment close to tangencies. A hint to this is already contained in the 
classical analogue representing a Poincare map from the boundary V to itself, which is 
discontinuous when the circle is a part of V. The same system can be described by a 
mapping which is when the circle is excluded from T and the mapping function 
includes the reflection the circle. 



2.1. The circle Green function 



In this section we consider a circle of radius R centered at the origin, with Dirichlet 
boundary conditions (^(r) = 0) on its boundary. In |Appendix A| we indicate how the 
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Figure 2. The different regions of r' for a given r. The penumbra occupies a small 
region on both sides of the geometrical shadow line (dotted) . 

results of this section modify for Neumann (dnip(r) = 0) or mixed boundary conditions 
(«V(r)+9ftV(r) = 0). 

The Green function of the circle satisfies (§) for any r, r' outside the circle, outgoing 
boundary conditions at infinity, and Dirichlet boundary conditions on the circle. The 
exact expression for the Green function is (in polar coordinates) 



G M =1 o E [Hr(kr < )+S l (kR)H+(kr < )]H+(kr > )e llAe 



(11) 



l=— oo 



where r> (r<) is the larger (smaller) of r and r', and A6 = 6 — 6'. The elements of the 
scattering matrix of the circle which is diagonal in the angular momentum representation 
are given by 

H-(kR) 



Si(kR) 



(12) 



H+(kR) 

We continue the discussion assuming r, r' > [1 + (kR)- 2 / 3 ]R, so that the Debye 
approximation may be used for H^(kr) and H^{kr') if / < kR. In addition we assume, 
without loss of generality, that < A6 < ir. 

The Poisson summation rule is used to express the Green function as 



G(r,r')= G (m) (r,r'), 



(13) 



m=— oo 
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Figure 3. The geometry of trajectories in the lit region (in the upper part) and of 
creeping trajectories in the deep shadow region (in the lower part). 

where 

■ r+oo 

G (m) (r, r') = - dl [#f (Jfer<) + S l {kR)H+{kr < )] H+(k r> )e UAe+27riml . (14) 

We discuss separately three regions for r and r' (see figure Q). The first is the lit region, 
in which the usual semiclassical results hold. It is obtained from the m = term in 
(|T3|), while the m/0 terms always describe the contributions of creeping waves. The 
second is the deep shadow region, in which only creeping waves contribute and the third 
is the penumbra, in which the leading order contribution, the m = term, comes from 
nearly tangent orbits. The lit and the deep shadow regions lead to known contributions 
of periodic orbits and the corresponding derivations will only briefly be described for 
the sake of completeness, and for understanding where the employed approximations 
fail when approaching the transition region. 

2.1.1. The lit region In the lit region the integral G^°'(r, r') gives the usual 
semiclassical result. For the scattering matrix element Si(kR) the Debye approximation 
is used for I < kR, and for I > kR it is approximated by 1. The integral is evaluated 
in stationary phase approximation. There are two saddle points, which relate to the 
two classical trajectories from r to r' (see figure ^). One is direct, and the other 
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reflects specularly from the circle. The Green function is then a sum of two terms. 
The contribution of the direct path is given by 

1 / 2i \ 5 

G d (r,r') = - [ — ) e ifcL , (15) 



4 \irkL, 

where L — I r — r'l , and the reflected path yields 



G r (r,rO = -i 



2i 



2 e ifc(L r +L r ,) Qg\ 

_7rfc(L r + L r / + 2L r L r //a)J ' 

where L r (L r /) is the distance from r (r') to the point of reflection such that the length 
of the reflected trajectory is L r + L T >. The impact parameter of the reflected path is b T , 
and a = \J R 2 — fof. 

The semiclassical result fails when the reflected trajectory becomes nearly tangent 
to the circle. To be precise, the semiclassical approximation holds provided that 
b T < [1 — (kR)~ 2 / 3 ]R. In other words the angle of incidence 9 should be smaller than 
tt/2 - (kR)~ 1/3 . This condition is necessary for the semiclassical approximation of the 
scattering matrix Si(kR) to hold at the saddle point, and thus defines the borderline 
between the illuminated region and the penumbra. 

2.1.2. The deep shadow region In the deep shadow region the contour of the integral 



(14) for m = may be closed in the upper half plane. The integral is then calculated by 
summing the contributions of the poles of the integrand which are just the poles of the 
S matrix. The first few poles are close to kR (figure §) and give a contribution which 
is evaluated using the transition region approximation for the Hankel functions. The 
position of these poles is given by 

l n « kR + e^ 3 x n (¥p) 3 , (17) 



where — x n are the zeros of the Airy function Ai(x) [23]. The residues of the S matrix 
at these poles are given by 

i 



c 



(kR 



27rAi'(-a; r) ) 2 V 2 



3 



The result is interpreted as the contribution of a creeping wave 0. It starts from r 
along a path which is tangent to the circle and creeps along the circle in one of the 
available creeping modes n until it leaves to r' along a path which is again tangent to 
the circle (see figure |]). The contribution of a creeping wave to the Green function is 




Figure 4. The complex I plane. The squares represent the poles of the S matrix. The 
bold line shows the contour from ctioo, going through kR and around the poles of the 
S matrix to 0200. The contour does not cross the two lines defined by (21)-(^2|), with 

J] -> ±71-/2. 



where L r (L r /) is the distance from r (r') to the corresponding point of tangency, and 
7^°) is the creeping angle. The length of the creeping trajectory is L r + L T > + j^R. The 
diffraction coefficient D\ is given for each creeping mode by 

2* = e*/<(^)% B . (20) 

The angular momentum along the creeping orbit in a creeping mode n is l n . Its real part 
is slightly larger than kR. The angular momentum also has a positive imaginary part, 
which is the result of the continuous decrease of amplitude as waves leave the circle in 
a tangent direction all along the way. 

The same calculation holds for all m 7^ terms of (|HD , for any location of r and r'. 
For m > 0, G^^r, r') gives the contribution of creeping waves going around the circle 
an additional m full times in the clockwise direction, such that 7 <m - 1 = 7^ + 2iim. For 
m < 0, the integration contour should be closed from below, and the result describes the 
contribution of creeping waves going in the anti-clockwise direction \m\ — 1 full times. 

The Airy approximation fllTf ), (|I5D fails when inserted into (|T9| ) if the creeping angle 
does not obey 7 ^> {kR)~ 1 ^. This condition guarantees that the contribution of the 
poles rapidly decays with n and defines the borderline between the deep shadow region 
and the penumbra. If it is not satisfied, the sum in fl!LTS| ) has to be extended to poles 
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which cannot be described by (|TT| ) and ( |T8D 



2.1.3. The penumbra region In the penumbra region, between the lit and the deep 
shadow regions, the above expressions fail. In this section we obtain expressions valid in 
the penumbra, following the methods used by Nussenzveig |22[ for studying the problem 
of scattering off a sphere. The contour of the integral (14) for m = is first deformed 
in the upper half complex plane (see figure [|), so that it goes from aioo, through kR 
and around the poles of the S matrix to <x 2 oo. The limits a\oo and (T200 are defined by 
demanding that the contour does not cross the lines asymptotically defined by I = a\l\, 
where 



a 



exp 



and 



77 = e In 



7T 

M 2 



21 



(21) 



(22) 



e/cr< 

The limits I — > aoo and e — > are taken simultaneously in such a way that 77 — ► ±7r/2. 
The deformation of the contour allows a separate treatment of the two parts of the 
integral fll4j) . For the first part we immediately obtain 
i 



dlH-ikr^H+ikr^e 



IA6 



0. 



(23) 



(TlOO 



as the contour of the integral may be closed from above and the integrand has no poles. 
The second part of the integral is split into two parts, which we call the direct and 
glancing parts of the Green function for reasons which will become clear. The result is 

i 



G d (r,r') 



(T 2 CO 



dlH+(kr)H+(kr')e 



iiAe 



(24) 



and 



r, r 



CTlOO 

i 



dlSi{kR)H+{kr)H+(kr')e 



iiAe 



(T20O 



dl [Si(kR) - l}H+(kr)H+(kr')e 



iiAe 



(25) 



Each of these terms will now be treated separately. 

In the direct part (24) the Hankel functions are replaced by their Debye 
approximation, and the integrand is evaluated around its saddle point. The only 
difference in the contribution of the direct term from that of the lit region expression is 
that now the limit of the integration must be accounted for, as it is close to the saddle 
point. The result is just the lit region expression (|15"D multiplied by a simple factor. It 
is given by 



G d (r,r') 



F(oo)-F(u) 



2i 



2i 

TtkL 



AkL 



(26) 
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Figure 5. The geometrical setup of trajectories in the penumbra. In the upper part 
r and r' are in the classically illuminated region (6<j > R), and in the lower part they 
are in the classically shadowed region. 

where F(x) — C(x) + iS(x) is the Fresnel integral, and 
/ kJ \ 5 

*= —7 (R-h). (27) 

The impact parameter is again denoted by b^, and we have z = v/r 2 — b\, z' = a/ r' 2 — b 2 d 
(see figure |5]). The Fresnel factor is in general a complex number. It equals | for exact 
tangency. Close to the border of the illuminated region it approaches 1, and it tends 
to at the border of the deep shadow region. Although this is the correct limiting 
behaviour, ( p6|) is restricted to the penumbra and does not necessarily represent the 
correct interpolation between the illuminated and the deep shadow regions. 

In the glancing part (25) the main contribution of both integrands comes from the 
vicinity of kR. The integrals are evaluated using the following approximations: The S 
matrix is approximated by the transition region expression 



where 



kR) 

The rest of the integrand is taken at kR, using the Debye approximation for the Hankel 
functions. The result of all these steps is that the only I dependence in the integrand is 



x=[-£ 5 Y(l-kR). (29) 
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in the argument of the transition region approximations. By a change of the integration 
variable the problem reduces to finding the value of the constant 

_ e"/ 3 f° Ai(xe- 2 "/ 3 ) J_ r°° Ai(ar) 
c ~ 2^ i /_ oo dx Ai(xe 2 -/3) + 2V3 _/ d;E Ai(xe 27ri /3) " ^ 30 ^ 

The first term is the complex conjugate of the second, as was shown by Rubinow and 



Wu , and after a numerical integration over the second term, we have 

0.996193019928. 
Finally, the glancing term of the Green function is given by 

GJr, r') = - — ( fc ^_ e ifc(^+L r ,+fl 7 (°))+^/ 3 ; (31) 
Air k\jL r L r i 

where L T (L T i) is the distance from r (V), along a line which is tangent to the circle, to 
the corresponding point of tangency, and 

ry(°) = A9 — arccos — arccos ■ (32) 

When bd < R, the values for L T , L T i and 7^°- ) are equivalent to the ones in the creeping 
case. When bd becomes larger than R (see figure |]), the points of tangency cross and 7^°- ) 
becomes negative. Note that in contrast to fll9D , equation ( pTj) contains no exponential 
damping. 

2.2. Diffraction corrections in the trace formula 

The contribution of periodic orbits to the density of states of the billiard may now be 
calculated. As was explained above, the contribution of orbits which bounce times on 
the exterior boundary V is found from the multiple reflection expansion (|T0|) using the 
circle Green function and performing the integrals over the billiard boundary excluding 
the circle in saddle point approximation. For each circle Green function the appropriate 
expression is used, depending on the positions of Yi and r i+1 . In the lit region we have 
G = G*d + Gr (1T5HT61), in the deep shadow region G = G c ([0J), and in the penumbra 



G = Gd + G g (|26| , |31|) . In the deep shadow region the deviation from ikL in the exponent 
in (|I~9"D is taken as part of the prefactor, as it is proportional to fc 1//3 and can be considered 
a slowly varying function for large k. Then in all cases the expressions for the Green 
function are of the form Ae lkL suitable for the saddle point approximation. After the 
saddle point integration a purely classical path yields the Gutzwiller contribution of 
an ordinary unstable isolated periodic orbit corrected with the product of the Fresnel 
factors for all segments traversing the penumbra. An important difference for creeping 
orbits is, that the action entering the Green function G(r, r') can be represented as a 
sum of two terms, one depending exclusively on r and the other on r'. This is obvious 
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from (0), (Pip and (^). As a consequence, the matrix of second derivatives entering 
the saddle point approximation to fllCf ) decomposes into K^> independent blocks for an 
orbit (P) with > 1 creeping segments. Finally one obtains an expression which is 
formally similar to the one derived by Vattay et al. for the contribution of a creeping 
orbit in the deep shadow region 

1 A K(P) 

d <£*\k) = —i-lm TT Gf cl) Gf iS) . (33) 
m dk xx j j 

The semiclassical Green function along the classical segment j connecting two 

points of tangency of the path to the circle is given by 

i 

~(scl) _ , 1 / 2l I 2 ihLj+iKH/2 (OA) 

where Lj represents the path length, /ij the number of conjugate points and fj the 
product of all the Fresnel factors along the segment. The monodromy matrix is defined 
such that Myi = dx' ± /d(j)" with and x± denoting the direction of the path and the 
coordinate normal to it, respectively. The Green function along a creeping segment is 
given by 



47r J2 n r n^ lln 7j deep shadow 

-2c(kRy/ 3 e ikR ^- i7r / e penumbra 



Gf lS) = { (35) 



and r in (|33|) counts the number of repetitions for orbits which are multiple traversals 
of a primitive orbit. 

Equation (|33|) gives the diffractive contributions in terms of geometrical information 
from creeping orbits. However, if one is interested in the contribution of a purely 
classical orbit which includes a very shallow reflection from the circle and must therefore 
be treated in the penumbra approximation, it is also possible to give a correction to 
the standard semiclassical expression without explicitly calculating the corresponding 
creeping orbit. For this purpose, the length of the creeping orbit is approximated by 
the length of the classical orbit since both orbits approach each other as the reflection 
becomes closer to tangency. The monodromy matrix along the classical segment of the 
creeping orbit is replaced by the monodromy matrix along the classical orbit excluding 
the almost tangent reflection which is incorrectly described in the standard theory. In 
terms of this reduced monodromy matrix M and the angle of reflection 9 the Gutzwiller 
amplitude for the contribution of the orbit is 



7T 



TrM - 2 + 



2M 12 


-1/2 l 


.Rcos 9 


1/2 


Rcos9 


7T 


2M 12 





(36) 
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where we have assumed cos^ <C 1. Comparing this to (^) we obtain the factor 

-=L={kR)-VW" , (37) 
V7rcosy 

by which the standard result for the contribution of the orbit is enhanced due to 
diffraction. 

Finally we would like to point out the different dependence on k of the contributions 
from various periodic orbits. The contribution of an isolated periodic orbit, with all 
segments in the lit region, is 0(k°), as is expected in the semiclassical approximation 
disregarding diffraction. For each segment in the deep shadow region, the contribution is 
multiplied by 0[exp(— CA; 1 / 3 )], where C depends on the creeping angle corresponding to 
this segment. For each glancing segment in the penumbra a factor 0(k^ 1 ^ 6 ) is obtained. 
It is worth noting that, unless the orbit segment is precisely tangent, it will eventually 
fall, for high enough k, either in the lit or in the deep shadow category. 



3. Penumbra diffraction in the spectrum of the annular sector 

The simplest billiard with a domain which is exterior to a circle, is the annulus defined 
by two concentric circles with radii R\ and R 2 - The annulus billiard is integrable, 
and its periodic orbits form one parameter families (or manifolds). The ratio R2/R1 
may be chosen such that two families of primitive periodic orbits, and their repetitions, 
will traverse the penumbra of the inner circle (one with direct segments and one with 
glancing segments). As a first numerical test of the expressions for diffractive periodic 
orbits, we consider the contribution of these orbits to the spectral density of the annulus. 
The annulus has the advantage that since it is integrable, its periodic orbits are easy to 
find, and its eigenvalues may be calculated with a high accuracy and for high energies, 
reaching the limit (kRi) 1 ^ 3 3> 1. As will be discussed, however, diffractive periodic 
orbits which bounce more than once on the exterior circle are not easily described by 
the theory presented in this paper because of a problem which is specific to the annular 
geometry. For this reason we will consider only the shortest diffractive periodic orbit of 
the desymmetrized annulus which bounces just once on the outer circle. The example 
of the annulus is thus limited, and in particular it does not allow to check whether the 
expressions for the Green function of the circle give the correct results after a saddle 
point integration. 

In this section we consider a desymmetrized annulus, the annular sector of angle 
a = it /4, and with R 2 /R\ ~ y/2 (see figure @(a)). In the annulus with this radii ratio, 
the square orbits of the outer circle are nearly tangent to the inner one. The orbits 
of the annular sector are obtained from those of the full annulus by desymmetrization. 
However, as a result of the symmetry of the square orbits, their desymmetrizations 
repeat the corresponding primitive orbits of the annular sector four times. There are 
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Figure 6. (a) The annular sector, with R 2 = V2Ri and a — 7r/4. Three T-orbits 
(which in this case are exactly tangent) are shown in dashed lines, (b) The absolute 
value of the length spectrum, with R2 = \pl and fco = 1300. Peak T corresponds to 
the T-orbits. Peaks D„ correspond to n repetitions of the orbits along the diameter. 
To peak P contribute orbits which do not reach the inner arc, and bounce n > 2 times 
off the outer one. 

two types of primitive diffractive orbits in the annular sector, the direct ones, arising 
from square orbits of the full billiard, and glancing orbits. They are referred to as T- 
orbits and bounce once off the external arc. For R 2 /Ri = V2 they merge into a single 
tangent orbit. 

The eigenstates of the annular sector are given by 

Vv(r, 9) = [Yiik^Jiik^r) - J^R^Y^r)] sm(W) , (38) 

for I — 4, 8, 12, . . . and n — 0, 1, . . . where k^ n are solutions of the secular equation 

Y l (kR 1 )J l (kR 2 ) - J l (kR l )Y l (kR 2 ) = . (39) 

We calculated the eigenvalues for R x = 1 and R 2 = \pl with k < 2800 (244 397 
eigenvalues) to an accuracy of 10~ 10 , and for several values of R 2 ~ \/2 in the region 
1050 < k < 1550 to the same accuracy. In order to extract the contribution of the T- 
orbits to the density of states we consider the truncated length spectrum (§) described 
in the introduction. The weight w(k) is taken in all calculations of this section to be a 
Gaussian with a center fco and a variance a 2 = 50 2 (see figure |6](b)). The peaks in |-D(x)| 
are Gaussians with a variance 1/cr 2 centered at lengths of periodic orbits. Our choice of 
the variance assures that the peak corresponding to the T-orbits is well separated from 
all other peaks. Thus we can study the contribution of T-orbits without considering 
any other orbit. 

To find the semiclassical contribution of the T-orbits in the annular sector, we 
consider a desymmetrized form of the circle Green function, which satisfies Dirichlet 
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boundary conditions on the straight lines of the annular sector. It is given by 

G(x, y, x\ y') = G(x, y, x', y') - G(x, y, x\ -y') - G(x, y, y', x') 

+G(x, y, y', -x') + G(x, y, -y', x') - G(x, y, -y', -x') 

-G(x, y, -x', y') + G(x, y, -x\ -y') . (40) 

The contribution of orbits which bounce once off the outer arc V to the density of states 
is found using ( |T0D as 



Im 



1 A. 

it dk 



ds 



dG 



(41) 



Different terms of the desymmetrized Green function (40) contribute to different orbits. 
The fourth and fifth terms are those which contribute to the T-orbit manifold. Their 
contribution is the same, as they correspond to the same orbit traversed in opposite 
directions. The integrand is independent of the integration variable and the contribution 
of the T-orbit manifold is therefore given by 



Im 



Rid_ 
2 dk 



d_ 

Or 



G(r,0; 0,r) 



(42) 



r=R 2 



For the Green function we use the expressions given in section |2.1| . It is then either of 
the form Ae lkL or a sum of two such terms. To leading order, which our calculations are 
limited to, both derivatives in (|42| ) act on the term e lhL only. 

In figure |7|(a) we present the results for the exactly tangent periodic orbit (i? 2 = \/2). 
The success of the penumbra approximation in this case is evident. There are two 
contributions, which are 0(k 1 ^ 2 ) for the direct path, and 0(fc 1//3 ) for the tangent path. 
The remaining error is very small, and decreases like fc _1//3 (the measured slope of the 
error line in the graph is -0.317), suggesting that the next order contribution comes from 
the glancing part of the Green function. This could be expected since the penumbra 
approximation which was used for the scattering matrix (|28| ) is accurate to Q(k~ 2 ^ 3 ). 

In figure f7|(b) we present the results for the peak from the T-orbits as a function 
of the external radius of the annular sector while /c is fixed. Starting from R 2 > R\, 
the semiclassical approximation is seen to break down as the orbit approaches tangency. 
On the other side the creeping approximation breaks down when the creeping orbit 
approaches tangency. The penumbra approximation is best at exact tangency. On both 
sides of the penumbra there is a region where none of the approximations is successful. 
Some form of uniform approximation would be needed to cover all regions well. Note 
that all the values of fed — R\ presented are inside the penumbra (6 r = 1 — k 2 ^ 3 for 
&d — R\ ~ 0.15, see section |2.1.1| ). The criterion v = 1 yields that the semiclassical 
expression for the contribution of the direct periodic orbit is valid when b^ — Ri > 0.035, 
and the error of this approximation can be seen to grow when bd — R\ becomes smaller 
than this value. 
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Figure 7. (a) The maximum of the T-orbits peak in |Z) w (x)| (full circles) as a function 
of ko for i?2 = V2 and R\ = 1 (exact tangency). The maximal error of the penumbra 
expression for this peak, including both direct and glancing contributions is given by 
the empty circles. Note the logarithmic scale, (b) The maximum of the T-orbits peak 
in |D w (x)| (full dots) for different values of i?2, with fco = 1300 and R% = 1. The 
maximal error of the semiclassical approximation for this peak is given by the crosses, 
for the creeping approximation by the squares and for the penumbra approximation 
by the empty circles. The whole region of bd — R\ is well inside the penumbra. In both 
graphs the lines are for guiding the eyes only. 



At the beginning of this section it was stated that there is a problem, unique to the 
annulus and the desymmetrizations thereof, in accounting for orbits which bounce more 
than once on the exterior circle. The length of an orbit with two consecutive creeping 
parts in the annulus is unchanged if the point on the exterior circle between these parts 
is allowed to vary (as long as both parts remain creeping). Thus, an orbit bouncing 
n times on the exterior circle and creeping in all its segments is in fact part of a n 
parameter family of orbits. The same situation occurs in the penumbra region for the 
tangent orbits. As this problem is unique to the annulus, we did not pursue it further. 



4. Diffraction effects in the quantized Sinai billiard 

In this section we study the Sinai billiard, i.e. an example for a chaotic billiard exterior to 
a circle. We demonstrate the importance of corrections to the standard Gutzwiller result 
due to the diffraction on the concave part of the boundary and check quantitatively the 
different approximations for the circle Green function derived in section |2.1|. Moreover 



we discuss also examples for orbits, where none of the described approximations accounts 
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for the diffraction corrections in a satisfactory way. 

The classical and creeping periodic orbits were calculated using the minimum and 
the unique coding principles |26|, PTfl . The quantum data on the Sinai billiard have 
been obtained using the scattering approach to quantization || . The scattering system 
we consider here has been described in the introduction (see figure [I]) and a detailed 
discussion of the numerical evaluation of the corresponding scattering matrix S(k) can 



be found in flTI 



However, the scattering approach to quantization does not only provide a framework 
for an efficient quantization of the Sinai billiard, it allows also to formulate the 
semiclassical theory in a way which is particularly well suited to study higher order 
corrections to the standard results. In order to demonstrate this we will study both the 
semiclassical density of states for the Sinai billiard, and the semiclassical approximation 
to the corresponding S-matrix. 



4-1. Analysis of the spectral density 

We begin with the energy spectrum of the quarter Sinai billiard with a = 1 and R = 0.5. 
The results are analyzed using the length spectrum (Q). Semiclassically, every periodic 
orbit (manifold) contributes to D w (x) in a small vicinity of its length, and this allows 
us to pinpoint individual contributions. The selected spectral interval < k < 300 
contains 5667 levels. The weight function w(k) was taken as a Gaussian centered around 
k Q = 150 whose width is a = 40. In the figures we show \D w (x)\ which is sensitive to 
both amplitude and phase deviations. 



4-1.1. Exactly tangent orbits In figure ^ one observes clear deviations between the 
quantum (exact) and semiclassical length spectra localized near the bouncing ball 
manifold at x = 2 and its double traversal at x — 4. The semiclassical expression 
contains the leading contributions from the bouncing ball families ||, the unstable 
isolated periodic orbits and the edge orbit Thus the deviations are mainly due to 
penumbra diffraction. As for the annulus, we use an appropriately desymmetrized circle 
Green function 

G[(x, y), (x , y')\ = G[(x, y), {x\ y')\ + G[(x, y), (-a/, y')\ 

+G[(x, y), (x', -y')\ + G[(x, y), (-x', -y')\ , (43) 

in order to calculate the penumbra corrections. In the multiple reflection expansion ([10|) 
we concentrate on the terms that give the bouncing ball contributions. For the shortest 
family and a single traversal we consider 

2 H f f\ dG[(x, yi ),(x,y 2 );k} 



dy 



y\=(-yi)=a, 
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Figure 8. Penumbra corrections of the length spectrum for the case of exactly 
tangent orbits. The orbits considered are drawn below the frames. Solid lines show 
the quantum (exact) length spectrum, and dotted lines show the difference between 
the quantum length spectrum and the semiclassical approximation supplemented by 
the results of [0 for the bouncing ball and edge contributions. Vertical bars indicate 
locations of unstable periodic orbits, daggers indicate bouncing ball families. Note 
the logarithmic scale, (a) Shortest tangent orbit (x = 2). Dashed line - glancing 
contribution also included, (b) Double traversal of the orbit considered in (a), x — 4. 
Dashed line includes 3 penumbra contributions (see text). 



+ 



dy 



dG[(xx,y), {x 2 ,y);k] 



(44) 



xi = (—x 2 )=a 



Due to the x «-> y symmetry, the two terms are equal, and it is enough to consider 
only one of them. We substitute for G its leading term approximation in the penumbra, 
G ~ Gd + G g and using equation (^) and the results of |Appendix B| , we get 



dbM(*) = <i(*) + dgw(*) 
db b l (Jfe) = — / cos ( 2ka - 
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(45) 
(46) 

(47) 
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The first term is the semiclassical contribution due to the bouncing ball family 0. 
The second term is a genuine diffractive contribution, that can be attributed to the 
exactly tangent orbit at the closure of the bouncing ball manifold near the circle. The 
contribution of this exactly tangent orbit is 0(k^ 1 ^ 6 ) which is slightly smaller than 0(k°) 
for unstable periodic orbits. 

In figure §(a) we present a portion of the length spectrum near x = 2. The 
differences between the quantum and the semiclassical predictions (including bouncing 
ball contributions) are about 10 percent. If we supplement the semiclassical expression 
with the tangent contribution in (47), the deviation reduces to about 1 percent. This 
clearly assesses the penumbra theory for exactly tangent orbits. 

A more complicated situation arises for the double repetition of the above bouncing 
ball family. The semiclassical contribution is given by an integral over a multiplication 
of two circle Green function, each decomposed into Gd + G g , which results in three 
terms: 



4 c* d f a j_ ^ dG[Pi,?/i), (x 2 , y 2 )\ dG[{x 2) y 2 ), {x x ,y x )\ 



4 d r 

4b,2(&) = / dx x dx 2 

71 dfc Jo 



dyi dy 2 
= 4Uk) + dil 2 (k) + dll 2 (k). (48) 

The terms are interpreted as "direct-direct" , "direct-glancing" and "glancing-glancing" 
contributions, with obvious notation. Straightforward calculations (see also [Appendix 
0) give 

d*i (jfc) = v^(a -fl) al/2fcl/2 cog r Aka _ ?r\ _ a cQs{4ka) (49) 

7T ' V 4/ TT 2 

4l 2 (k) = - ^a^RWk-W cos (ika + ^) (50) 

2 

rfffi Jk) = ^-R 2 /^- 1 / 3 cos ( Aka + - ) . (51) 

' 7T V 6/ 

The results for d[\ h v d^ 2 and d^ 2 agree with the general expression (|SB| ) and can be 
interpreted as the contributions from isolated orbits which are exactly tangent to the 
circle. The term d^ 2 contains the semiclassical contribution of the bouncing ball family 
and an interesting correction which is of the same order as an unstable periodic orbit. 
This correction comes from the non-zero average of the squared Fresnel factor and thus 
is of diffractive origin. There is no unstable periodic orbit that gives this contribution. 
These predictions are fully verified against the numerical data shown in figure |§(b). 
Indeed, the penumbra corrections reduce the deviations very significantly near x = 4. 



4-1.2. Almost tangent orbits We now turn to the investigation of almost tangencies, 
that occur in generic isolated and unstable periodic orbits. In figure | we indicated a 
few periodic orbits for which there are significant deviations due to penumbra effects. 
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Figure 9. Penumbra corrections of the length spectrum for the case of almost tangent 
and ghost orbits. The orbits considered are drawn below the frames. Solid lines show 
the quantum (exact) length spectrum, and dotted lines show the difference between the 
quantum length spectrum and the standard semiclassical approximation. Vertical bars 
indicate locations of unstable periodic orbits, (a) Pair of almost tangent periodic orbits 
ati« 5.10. Dashed line - direct term included, sparse dotes - glancing contribution 
also included, (b) Pair of classically forbidden periodic orbits at x ~ 5.24. Their 
location is indicated by a double dagger. Notation is as in (c). 



We choose to concentrate on the pair of periodic orbits of lengths 5.10845 (direct) and 
5.10908 (glancing) which are plotted in figure |9](a). The orbits are geometrically similar, 
except that the direct orbit has a segment that just misses tangency with the circle, while 
for the glancing orbit the corresponding segment reflects from the circle in a very shallow 
angle. This angle of reflection is about 1°, which is well inside the penumbra and fully 
justifies the implementation of the corrections. To calculate the corrected contributions, 
we replaced the semiclassical Green functions with their penumbra counterparts Gd and 
G g (see ( p^) , (pip). For the direct orbit, the only change was a multiplication of the 
standard semiclassical contribution by by a Fresnel factor, whose value was 

F(oo) - F(v* -0.35) ^ Q 71e - .23i 



for k = ko = 150. Including this correction reduces the deviations significantly, as can 
be clearly seen from figure |5](a). To account for the glancing corrections we use the 
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approximation (|3T|). While the difference in the lengths of the classical orbit and the 
corresponding creeping orbit (with negative creeping angle) is very small (« 10" 6 ), the 
prefactor significantly grows by a factor of ~ 4.7. The contribution of the glancing 
orbits further reduces the deviation by a factor of 2, as seen in the figure. 



4-1.3. Ghost orbits One of the most interesting applications of penumbra corrections 
is for ghost and creeping orbits, which are classically forbidden. Ghost orbits which are 
almost tangent in the shadowed part of the penumbra, are expected to give appreciable 
contributions, comparable to standard semiclassical contributions of real periodic orbits 
with similar lengths. To find such ghost orbits in the Sinai billiard, one needs therefore to 
look at periodic orbits that are pruned at a radius slightly smaller than R = 0.5 which 
we use for the quantum results. Indeed, we observe a pair of geometrically similar 
periodic orbits that coalesce and prune at R ~ 0.48. After enlarging the radius back to 
R = 0.5, we get a pair of "direct-shadowed" and "glancing- shadowed" penumbra orbits 
(see figure @(b)). The lengths of the orbits are ~ 5.2409 and 5.2413, respectively. The 
creeping angle of the glancing orbits is ~ 1.4°, which is small enough to justify the 
penumbra approximation. The Fresnel parameter v for the direct orbit is ~ 0.25, that 
gives a multiplicative Fresnel factor « 0.39 exp(0.31 i) which indeed indicates almost 
tangency. The direct contribution is by a factor « 3 larger (and with opposite sign) 
than that of the glancing orbit, and thus we should expect to see a noticeable peak in 
the length spectrum. Our expectations are fulfilled, as can be seen in figure ||(b). We 
can identify a peak in the quantum length spectrum near x = 5.24 with large deviations 
between the quantum and the standard semiclassical results. They correspond to the 
ghost orbits, and if we include ghost contributions, the deviations significantly decrease, 
which indicates the success of the theory. We tried to naivly implement the geometrical 
theory of diffraction as in [|7|, which takes into account only the creeping orbit. Summing 
over many creeping modes to get a convergent answer, we obtained large deviations from 
the quantum results in the length spectrum, as expected due to the small creeping angle. 



4- 2. The semiclassical S-matrix 

In this section we consider the semiclassical approximation to the S-matrix involved 
in the scattering quantization. As discussed in the introduction, all the spectral 
information on the billiard is contained in the total phase and all the traces of the 

5- matrix. We would like to show how this information may be extracted and used to 
study very fine details of the spectrum which would otherwise not be accessible. The 
fact that the S-matrix is a continuous function of k rather than a sum of delta peaks 
necessitates an analysis which is slightly different from that presented in the sections [3] 



and [O. The central spectral quantity which we consider here is the number counting 
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function N(k) = J k dk' d(k') rather than the density of states. According to (|J) it can 
be decomposed as 

oo 

N(k) = NW (k) + N(n) ( k ) ( 52 ) 

n=l 

with 

N( R \k) = ^Q(k) , 
2tt 

53 

(—1) n 

N {n) (k) = ± ^ImTrm) n . 

nn 

Each of the terms in this decomposition can be analyzed separately. We will concentrate 
here on the first two terms and . N( R ' will be referred to as resonance counting 



function for reasons explained in [28|. Beside the oscillating contributions from trapped 
periodic orbits it contains also a smooth part which is identical to the smooth part of 
the billiard spectrum up to the third term in the expansion of Weyl's law . This is not 
completely obvious, since it is known for the scattering from the outside of a billiard 
|20| that in general the above statement holds only for the area term. For the Sinai 



billiard we have clarified this point in the introduction by considering the limiting case 
of an empty square billiard where an analytic expression for the total phase is available. 

In order to pinpoint the contributions from individual periodic orbits we consider 
again length spectra which are now obtained by a fast Fourier transform using the 



Welch window and based on discrete points with a spacing Ak. Note that due 
to the discreteness of the transformation the contributions from long orbits appear at 
x = 

L (P) 

mod27r/Afc and may interfere with the contributions from the shorter orbits of 
interest. An example is provided in figure [10|, which displays the length spectra obtained 
from the oscillatory parts of N(k) (top), N^(k) (middle) and (bottom) using a 
logarithmic scale. From the number counting function we obtain a few well pronounced 
peaks at short lengths which are on top of a very large background containing the 
contributions from all long orbits. Although these long orbits are extremely unstable, 
the exponential proliferation of orbits ensures that the amplitude of their combined 
effect diminishes only slightly with the orbit length. 

The situation is different for the trace and the total phase of S, where the number 
of contributing orbits does not grow exponentially with the length. Nevertheless the 
long orbits become exponentially unstable and this is why we see almost no background 
in the lower two parts of figure So the observation of very fine structures in the 
length spectrum becomes possible, e.g. the tiny peak at x ~ 3 in the spectrum of the 
resonance counting function, which is due to diffraction effects as we shall explain in 
the sequel. It would hardly be possible to discover and study such peaks among all the 
leading order contributions from the unstable periodic orbits with a similar length. The 
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Figure 10. Comparison between the length spectrum of the number counting function 
N(k) (top), the term N^(k) related to the trace of the S-matrix (middle) and the 
resonance counting function N^ R \k) (bottom). The length spectra are obtained from 
the S-matrix for a = 1 and R = 0.5 in the interval k/n = 1, . . . , 281 (containing more 
than 49 200 eigenvalues) and with A fe = n/64. The smooth parts of N(k) and iV( R )(fc) 
have been subtracted according to Weyl's law. 



27 



u c 



e c 



r.12) 

^2 



(3) ,,(3) 



Figure 11. The figures show some trapped orbits of the opened Sinai billiard, bb 
denotes an example for a bouncing ball orbit and e is the edge orbit limiting the family. 
The other edge orbit u is the only isolated classical orbit contributing to the resonance 
counting function. The orbits denoted with c are diffractive. 
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Table 1. The table gives the necessary geometrical data for the computation of the 
contributions from some diffractive orbits to N^(k). L< cl ) and M[f } are the path 
length and monodromy matrix of the classical segment of the orbit, 7 is the creeping 
angle. 



resonance counting function is particularly well suited for a semiclassical analysis, since 
it contains contributions only from the very sparse set of trapped orbits, i.e. those orbits 
which never hit the section Tg of figure [l]. Some of these orbits - classical and diffractive 
- are displayed in figure [11. 

The most important of the oscillatory contributions is of order k 1 ! 2 and comes from 
the family of neutral bouncing ball orbits bb with length 2a and its multiple traversals. 
It is responsible for the large peaks at x = 2 and x = 4 in the length spectrum of the 
resonance counting function which is displayed in fig. [12] for a = 1 and R = 0.5. The 
dashed line is obtained by subtracting Weyl's law from the exact quantum result and 
contains therefore all oscillatory contributions to N^(k). 

Disregarding diffraction but including the edge orbit e running along the right 
billiard wall, the contribution of the bouncing ball family has been derived in and is 
also contained in equation (0). Beside this, the only standard semiclassical contribution 
to the total phase comes from the isolated unstable orbit u running along the left billiard 
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Figure 12. The length spectrum of N^^k) obtained from the interval 1 < ka/n < 
281 for a billiard with a = 1, R = 0.5 and Dirichlet boundary conditions on the circle. 
The dashed line corresponds to the full quantum result with the smooth part subtracted 
according to the generalized Weyl law, i. e. it comprises all oscillatory contributions 
to the resonance counting function. The thin solid line shows the deviations of the 
semiclassical approximation based on all classical orbits from the quantum data and for 
the thick solid line the diffractive orbits have also been included into the semiclassical 
approximation. The lengths of contributing orbits are marked with vertical bars. 



boundary. The semiclassical amplitude of such an edge orbit has been derived e.g. in 
]27J . For the thin solid lines in figure |T2| all standard semiclassical contributions have 
been subtracted. The remaining peaks in the solid line are now exclusively due to higher 
order corrections to Gutzwiller's result. At the length of the unstable orbit u (x = 1) this 
is sufficient to reduce the amplitude of the peak by more than two orders of magnitude. 
The subtraction of the leading order semiclassical expression is less successful at the 
bouncing ball lengths, since there we have very large diffraction corrections, which are 
explicitly given in (47) and (49)-(51). All these terms have been subtracted from the 
quantum mechanical data to obtain finally the thick solid line. Indeed, the magnitude 
of the peaks at the bouncing ball lengths is now also reduced considerably. 

The analysis can be further supported by considering the dependence of the 
magnitude of the peaks on k as displayed in figure [13] using a double logarithmic scale. 
The curves are obtained by restricting the Fourier transformation of the data to small 
intervals centered around a mean k given in units of n/a at the abscissa of the plots. 
The width of the intervals is chosen large enough to guarantee a sufficient resolution in 
the length spectrum. The approximate exponents a for the fc-dependence of the peaks 
in the semiclassical domain have been determined by a linear fit to the curves at high 
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Figure 13. The dependence of some peaks in the length spectrum of (k) on the 
wave number A; in a double logarithmic plot. The length spectrum was calculated in 
k- intervals of width 307r/a centered at the value which is given on the abscissa in units 
of n/a. As in figure [l2] the upper curve corresponds to the quantum result while the 
middle and the bottom curve show the deviations of the semiclassical approximation 
excluding and including diffraction, respectively, (a) peak at the length of the unstable 
isolated orbit u. (b) and (c) one and two traversals of a bouncing ball orbit. The k- 
dependence can be fitted by a power law with an exponent given next to the curves. 



values of k and are given next to the curves. In figure [15(a) we demonstrate in this 
way that the contribution from a standard periodic orbit is correctly described by the 
Gutzwiller formula up to corrections of the order k" 1 . At x = 2 (b) and x — 4 (c) the 
dashed line corresponding to the oscillating part of the resonance counting function is 
close to a fc 1//2 -behaviour which indicates that the peak is dominated by the bouncing 
ball contribution. The exponent for the peak at x = 2 in the solid line is close to 
— 1/6 which agrees with the penumbra contribution of the tangent orbit (47). When 
this expression is also subtracted from the data (fat line) we are left with a peak whose 
magnitude is much smaller and moreover falls off faster than A; -1 / 2 , i.e. the description 
of the tangent orbit is quantitatively correct and further corrections are of lower order. 
The situation is similar at x — 4. However the leading order correction to the bouncing 
ball result is now ~ k° and does not come from an isolated tangent orbit but from the 
diffraction correction in (49) to the family itself. When all diffraction effects have been 
subtracted (fat line) the height of the peak is reduced by an order of magnitude, but 
the peak does not fall off as fast as for x = 2. This apparent discrepancy is at least 
partially due to the finite available fc-range since it disappears gradually as the interval 
of computation is further enlarged (not displayed). We also recall that another error 
may be due to the explained "folding" of the peaks from long orbits on top of the length 
spectrum. 

Beside the discussed corrections to the standard semiclassical contributions, the 
solid line in figure [12| displays a number of additional peaks which are exclusively due 
to diffractive orbits. The open Sinai billiard supports four families of primitive creeping 
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Figure 14. The length spectrum of N^(k) computed from the interval 1 < ka/n < 
81 for a = 1, R = 0.5 and Neumann boundary conditions on the circle. The line types 
are chosen as in figure O. 



periodic orbits. The two shortest members from each of these families are displayed in 
figure [□]. These are orbits which are not due to a bifurcation of the type described in 
the introduction or in the last section. Rather they remain creeping even for very small 
radius. The creeping length of the orbit Cj° is exactly zero, and therefore this orbit has 
always to be treated within the penumbra approximation and results in (47). The other 
orbits are in the penumbra or deep shadow region depending on the value of k. Table [TJ 
contains the necessary geometrical information to evaluate the contributions from the 
shortest members of the four families. In figure ^ the lengths of the creeping orbits 
are indicated with vertical lines and the most important of them are denoted on top 
of the plots. It is clearly seen that each peak in the length spectrum corresponds to 
a particular periodic orbit. The parameters are such that the creeping orbits are well 
described by the deep shadow approximation throughout the whole k range. Indeed 
the magnitude of the peaks in the thin solid line is reduced when the contribution of 
these orbits is subtracted according to fl33|) (fat curve), although the reduction is not as 
striking as for the exactly tangent orbits. 

From the expressions derived in the Appendix A we expect that in the case of 
Neumann boundary conditions the creeping orbits have a larger amplitude than in the 
case of Dirichlet boundary conditions. Indeed the peaks at the corresponding lengths 
in figure [13] (Neumann b. c.) are more pronounced than in figure [12] (Dirichlet b. c.) 
and the success of the deep shadow approximation for these orbits is even more evident 
(note the different scaling of the ordinate according to the different magnitude of the 
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Figure 15. The length spectrum of N^(k) computed from the interval 1 < ka/n < 
81 for a = 1, R = 0.8 and Neumann boundary conditions on the circle. The line types 
are chosen as in figure |l2|. 



peaks). It is also possible to observe and correctly account for a peak at x = 7r/2 due to 
the circumference orbit C which never leaves the circle at all. As the radius of the circle 
grows, the lengths of the creeping orbits become closer to each other and the length 
spectrum is more complex (figure |T3]). Nevertheless we see that the magnitude of the 
peaks can be reduced considerably, when the semiclassical contributions are subtracted 
according to (|3"3"D. Unlike classical orbits, an arbitrary combination of primitive creeping 
orbits can be joined to form a new periodic orbit which also gives a contribution to the 
spectrum. Of particular importance are the combinations including the tangent orbit 
c$ , since then the creeping angle is relatively small. An example for a creeping orbit 
of this type is the peak denoted with + c± in figure [15. 

Now we would like to discuss the semiclassical approximation to Tr S in more detail. 
Again we will show that important corrections to the leading order semiclassical result 
are due to diffraction effects. In figure 16 we display the length spectrum of Tr S, which 
can be expressed semiclassically in terms of all periodic orbits of the Sinai billiard hitting 
the section Tg exactly one time. This includes the contributions from the two bouncing 
ball families present in the Sinai billiard with R = 0.5 which result in the two most 
prominent peaks in the spectrum at x = 2 and x = 2a/2. Most of the other large peaks 
can be related to isolated unstable periodic orbits as it can be seen from the vertical 
bars in the upper part of the figures which denote the lengths of all the classical orbits 
contributing to S. The remaining difference between the semiclassical approximation 
and the quantum data after the leading order contributions from all classical orbits 
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Figure 16. The length spectrum of Tr S(k) for a = 1 and R = 0.5 in the interval 
ka/ir — 1, . . . ,281 with = 7r/64. The full line shows the quantum data and the 
shaded areas represent the error of the leading order semiclassical approximation. The 
location of the periodic orbits contributing to Tr S is marked with vertical bars. 



have been included according to the standard Gutzwiller result is displayed using the 
grey shaded areas. We see that the semiclassical result accounts very well for some of 
the peaks at small lengths: all the peaks up to x = 4 are reduced by the subtraction 
of the semiclassical result by at least a factor of 10 (note the logarithmic scale). The 
result becomes increasingly worse as the length grows and from x ~ 7 the leading order 
semiclassical theory fails completely. 

In order to explain where the deviations come from we have further enlarged the 
region 4 < x < 5 in figure [17]. There are 12 isolated classical orbits in this interval 



which are displayed in figure 18. We observe large deviations between the semiclassical 
result and the quantum data for the orbits with L = 4.12,4.22,4.61 and 4.62. These 
orbits have in common that one of their reflections from the circle is very shallow, 
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Figure 17. The length spectrum of TrS(k) for a — 1 and R = 0.5 in the interval 
ka/ir — 1,...,281. The full line shows the quantum data and the shaded areas 
represent the error of the standard scmiclassical approximation. The location of the 
unstable periodic orbits contributing to Tr S is marked with the dotted vertical lines. 
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Figure 18. The unstable periodic orbits contributing to TrS'(fc) with lengths in the 
interval 4 < L < 5 corresponding to figure [l7]. Below the figures the length and the 
stability prefactor of the individual orbits are given. 
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Figure 19. Some of the orbits which give rise to the peaks at large x in the length 
spectrum of Tr S. Below the orbits the length and the dominating eigenvalue of the 
monodromy matrix are given. Due to diffraction the standard Gutzwiller formula fails 
to predict the contributions of these orbits correctly. 



i.e. the classical description excluding diffraction effects approaches the limits of its 
validity. However, unlike the orbits discussed in the context of the spectral density, 
where a considerable improvement of the semiclassical result could be achieved when 
the penumbra corrections are taken into account, the orbits of figure IT3 are not really 
well inside penumbra but rather at the border between the illuminated region and the 
penumbra. Consequently the results of section |2| do not apply. As already mentioned 
at the end of section |3|, a uniform geometric theory of diffraction would be needed in 
order to correct the contributions from such orbits. 

The peaks at 4.47 and 4.75 do not correspond to classical orbits but represent the 
contribution from creeping and "ghost" orbits. L = 2y/5a = 4.472 is the length of 
the first completely shaded bouncing ball family, which does not give a contribution to 
leading order. However, an important diffraction correction resulting in the observed 
peak in the length spectrum is due to those orbits in the family which traverse the 
penumbra. 

Finally we turn to the conspicuous peaks in the length spectrum of Tr S which 
are located in the vicinity of x — 8, 10, 12, . . . and can be explained by the orbits 
displayed in figure [B5[ Below the orbits the length and the dominant eigenvalue of the 
monodromy matrix are given. The orbits become very unstable as the length increases 
and the standard Gutzwiller expression fails completely to predict the amplitude of their 
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contribution, which is due to the almost tangent reflection from the circle. Although the 
penumbra approximation is much better and predicts at least the order of magnitude of 
the contributions, it is not capable of giving a satisfactory quantitative description and 
therefore not displayed. While for the shorter orbits the reflection from the circle is not 
yet inside the penumbra, the problem with the longer orbits is that more and more of 
the straight segments are close to the circle and need an additional diffraction correction 
which leads to an increasing error. Moreover additional corrections due to intermittency 
which were recently derived in [51 may be necessary to predict the amplitudes correctly 
since the long orbits are very close in phase space to the family bb of bouncing ball 
orbits transversal to the channel. 



5. Frequency of penumbra traversals 

The numerical examples of section HTT| illustrated the success of the theory derived in 
section |] to account for the significant penumbra corrections for particular periodic 
orbits. A natural question would be: how many periodic orbits should be corrected in 
this way? To answer this question we should consider two factors. First, the borderlines 
of the penumbra are k dependent and the fraction of phase space occupied by the 
penumbra is of order (kL)~ 2 ^ 3 , where L is a typical length of the billiard. This means 
that the penumbra shrinks to as k — > oo. In particular, we can also conclude that 
any given periodic orbit (that has no exactly tangent segments) will be outside the 
penumbra for k large enough. The second factor to consider is that in order to quantize 
the billiard up to a wavenumber k with a resolution of the mean level spacing one needs 
to consider periodic orbits up to the Heisenberg length Lh ~ kL 2 , or, equivalently up 
to number of bounces % ~ kL. This enhances the chance to visit a given area of phase 
space as k grows. To obtain the overall effect of these two contradicting trends, let us 
consider for each periodic orbit A = min,, \lj — kR\, where lj is the angular momentum 
of the jth segment of the orbit with respect to the circles center. The orbit traverses 
the penumbra at least once if A < (kR) 1 ^ 3 . If we assume ergodicity, then each segment 
of the orbit has an a-priori probability 

p « {kL)- 2 ' 3 (54) 

to traverse the penumbra. Assuming statistical independence of the segments, and 
homogeneous coverage of phase space by long periodic orbits, the probability that an 
orbit with n bounces avoids the penumbra is 

(1 - p) n « exp [-n{kL)- 2 ' 3 ] = exp ( ^ j . (55) 

Because of the exponential proliferation of periodic orbits, the overwhelming majority 
of periodic orbits satisfy n ~ % and thus n^ 3 < n. This means according to 
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Figure 20. The coarse-grained distribution of A/(kR). Note the logarithmic scale. 



that in the semiclassical limit most of the periodic orbits traverse the penumbra at 
least once, and for them the semiclassical approximation fails and should be corrected. 
To emphasize this point we rephrase our findings about the semiclassical limit in the 
following way: 

• Given a periodic orbit (which is not exactly tangent), its standard semiclassical 
contribution is recovered for k large enough. 

2 /3 

• For a given k, most periodic orbits which are shorter than m n H are described by 
the standard semiclassical approximation. 

• Although their number grows with k, their fraction out of the relevant periodic 
orbits becomes smaller, and the great majority of periodic orbits are affected by 
penumbra corrections. 

To verify these ideas we calculated A/(kR) for all periodic orbits of length up to 10 
of the quarter Sinai billiard. There were 20 150 primitive orbits, with total of 320 002 



segments. The coarse-grained distribution is shown in figure |20|. It is sharply peaked 
near the minimal value A = 0, which indicates that almost all periodic orbits include 
a nearly tangent chord, as predicted. To get a quantitative estimate, let us take the 
Heisenberg length to be the maximal periodic orbit length: Lh = 10. We need to 
estimate the relevant k. We use the definition of the Heisenberg time 



T H (E) = hd(E), 



(56) 
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which for billiards can be written as 

2nd(k) = L u (k). (57) 

If we use the leading order expression for d(k) for billiards 

Ak 

d(k) = — , (58) 
where A is the of the billiard we get 

k = ^ « 12.4 (59) 
in our case. The relative phase space area occupied by the penumbra is estimated as 

(60) 

where b max = y/2 is the maximal impact parameter in the billiard. The factor 2 is due 
to taking into account impact parameters which are both larger and smaller than R. 
The number of bounces is estimated by 

n R = — = — T ~ 15, (61) 
c ttA 

where c is the mean chord and 7 is the billiards perimeter. Due to the exponential 
proliferation of orbits, their majority will have a length close to Lh and thus hh chords. 
The probability of such orbits to avoid the penumbra is 

q = (1 -p)™ 11 w 0.03 (62) 

and consequently about 97% of the orbits are expected to traverse the penumbra. The 
penumbra borders in terms of A/(kR) are estimated as 

A 2 

— < (kR)s w 0.3 (63) 
kR 

which includes according to the numerical data 95.7% of the orbits, in good agreement 
with the theory. 



6. Conclusions 



In this paper we have considered the semiclassical quantization of billiard systems. Using 
asymptotic approximations to the circle Green function we have derived corrections 
to the standard Gutzwiller formula which account for the quantum diffraction from 
concave parts of the billiard boundary in the penumbra. This is the nearly tangent 
parameter region where neither the leading order semiclassical result nor the deep 
shadow approximation || [7|] are valid. There are two types of corrections: the first 
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can be expressed in terms of nearly tangent classical periodic orbits (including "ghost" 
orbits which cut straight through the concave billiard boundary). The contributions 
from these orbits differ from the Gutzwiller result by a prefactor of the order one, 
i.e. the correction can be as large as the standard semiclassical amplitude itself. The 
other type of correction can be expressed in terms of creeping orbits including those with 
negative creeping angles. Although these orbits are also obtained from the deep shadow 
approximation, their contribution is different in the penumbra: the amplitude does not 
decay exponentially in k but only as A; -1 / 6 for each creeping segment of the orbit. The 
appearance of this new type of orbits in the semiclassical quantization formulas raises 
questions about the structure of the set of all creeping orbits, e.g. how they can be 
computed using an extremum principle together with a code and if there is a one to one 
correspondence between the nearly tangent classical orbits and the creeping counterparts 
by which they have to be replaced in the case of penumbra diffraction. Although we 
have some preliminary results for the Sinai billiard, the answer to these questions is 
beyond the scope of the present paper. 

The derived corrections to the standard Gutzwiller result have been tested in both, 
the integrable annular billiard and the chaotic Sinai billiard. The success of the theory 
was obvious as long as the orbits under consideration were indeed well inside the 
penumbra region. 

In the case of the Sinai billiard we have used beside the length spectrum obtained 
from the set of billiard eigenvalues also an alternative method of spectral analysis which 
is based on the scattering approach to quantization. Here discrete length spectra are 
directly obtained from the total phase or the traces of the involved S-matrix which has 
the advantage that only limited subsets of periodic orbits contribute. The resulting 
sparse spectrum is particularly well suited to observe deviations from the standard 
Gutzwiller result and we could in this way check the quality of the penumbra and the 
deep shadow approximation to a high precision. 

The importance of penumbra diffraction corrections becomes obvious if one 
estimates the number of orbits which are affected by them. It turns out that the 
amplitudes of most of the orbits contributing to the spectrum up to some fixed value k 
must be corrected. In view of this fact it would certainly be desirable to further extend 
the results which we presented in this paper to the parameter regions, where none of 
the so far known expressions is sufficiently accurate. 
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Appendix A. The circle Green function with Neumann and mixed 
boundary conditions 

The different approximations to the circle Green function are found in a manner similar 
to section |2.1| when homogeneous boundary conditions are imposed on the circle. In 
this appendix we give the modifications to the expressions derived in section |2.1| for this 
case. Equations ( |i~3"l - |i4] ) are still the starting point for all calculations, and the boundary 
conditions affect only the scattering matrix Si(kR). 

For Neumann boundary conditions (daip(r) = 0) the scattering matrix is given by 

S l{ kR) = -%^±. (Al) 
V ; H+ (kR) V ; 

In the lit region, Gd(r, r') (|T5|), the contribution of the direct path, is unaffected, while 
the contribution of the reflected path G r (r, r') ([16]) is multiplied by —1. In the deep 
shadow region, the contribution of the creeping path G^\r, r') is still given by (|T^-|2Ti|), 
but with different poles l n and residues r n . The poles are given by (p~T|) , where in this 
case x n are the zeros of Ai'(x). The residues are given by 

e" i7r / 6 (kR" 

Tn ~ 27TX n Ai{-X n ) 2 \ Y 

In the penumbra the contribution of the direct path Gd(r, r') (|26|) is unaffected by the 
boundary conditions. The constant c in the contribution from the tangent path G g (r, r') 
(|3T|) is now given by 

1 f° , Ai'(xe- 2 ^/3) e ~2ni/3 ,00 A j/, ) 



(A2) 



2V3 y_ oo Ai'(xe 27ri /3) 2 V3 J Ai'(xe 2 ^ 3 ) ' 

The first term is again the complex conjugate of the second ||25|| , and by a numerical 
integration over the second term we have 

c w -0.864251443481 . (A4) 

We now consider the general case of mixed boundary conditions (Ki/j(r) + dhip(r) = 
0). Neumann boundary conditions are the case of k = 0, and Dirichlet boundary 
conditions are the limit k — > oo. The semiclassical quantization for mixed boundary 
conditions is considered in p0|. The scattering matrix is given by 



= KHr(kR)-kH r '(kR) 

lK ' KH+{kR)-kH+\kR)' 1 J 

In the lit region G&{r, r') (TTST) is again unaffected. The contribution of the reflected 
path, G T (r, r') ( P^) is multiplied by a phase e 1 *, where is given ]30| by 



2arctan I — cos 9 ) . (A6) 
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Note that for k — > oo (Dirichlet) = 0, and for k = (Neumann) <p = tt. The poles Z n 
for the contribution of the creeping path Gc (r, r') (|19| - p0| ) in the deep shadow region 
are still given by (|17]), where in this the solutions of 



Ai' 



a MT/3 



Ai(-x) 
The residues are given by 



k ( kR 
k 



-«7r/6 



2tt [Ai'( 



+ x n Ai( 



(A7) 



(A8) 



In the penumbra Gd(r, r') fl26|) is again unaffected. The constant c in G g (r, r') (|3l~l) must 
be replaced by the fc-dependent expression 



kR 



(A9) 



After some manipulations c(z) is given in a form convenient for numerical integration 
by 



c(z) 



2 l/3 



da; 



zAi(ar) + e-^ i/3 Ai'( 



^Ai(e- 2,ri / 3 x) + e 27ri / 3 Ai'(e- 27ri / 3 x) 



2V3 



dx 



zAi(z) + Ai'(x) 



(A10) 



zAi(e 27ri / 3 x) + e 27ri /3Ai'(e 27ri /3 x ) ' 

For some applications it is convenient to express mixed boundary conditions as 
6cosa"0(r) + smac^^r) = 0. The parameter a interpolates between Dirichlet (a = 0) 
and Neumann (a = vr/2) boundary conditions. In |30| the derivative 



d 

—d(k;b,a) 



(All) 



a=0 



is introduced as a tool for analyzing the spectrum of the Sinai billiard. It was conjectured 
there that the derivative of tangent contributions should semiclassically vanish. Using 
(A10) we find that 
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'kR' 




— c| 






cot a J 


da 




_ 2 _ 





3 27ri/3 



k 



a=0 



b {kR) 1 / 3 



(A12) 



which together with (47), (50) and (51) indeed gives a contribution 0(k 1//2 ) smaller 
than for a standard unstable periodic orbit. 



Appendix B. Calculation of direct contributions 

In this appendix we calculate the purely "direct" contributions to d^i(k) and d\^^{k) 
that appear in section [4. 1 . 1| . The purpose is to find genuine diffractive (penumbra) 
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contributions beyond the semiclassical ones 0. There is an inherent difficulty in this 
problem, since diffractive effects are localized near the circle, while the bouncing ball 
family is more global and covers a considerable volume of configuration space. In the 
lack on a uniform approximation for the direct contribution it is natural to consider 
a splitting of the integration region into "near" and "far" regions, using the appropriate 
expressions in each region. However, in our case it is not completely clear what is 
the correct splitting, and what is the correct transition between the penumbra and the 
illuminated and shadowed regions. Thus we choose to use the penumbra approximation 
of Ga for the whole integration region, and after performing the calculation to consider 
more closely the origin of corrections, if any. 

To get d bbl (k) we substitute in equation (44) the explicit form of Gd (equation 
([26])), and to leading order in k we have 




<iW = ^Ps<|e i ^f / F[2 X I— (/?-,-) \d.r} . (B13) 

7T2 

where we denoted F(x) = [-F(oo) — F(x)]/V2i. The integral is evaluated by using the 
indefinite integral of the Fresnel function 

h(t) = I* F(Odf = tF(t) - \\jl^ t2 . (B14) 

It can be simplified, if we consider the asymptotic approximation of the Fresnel function 
for \t\ > 1: 

F(t) = F(oo) sign(t) - — e^' 2 + 0(|tr 3 ). (B15) 

nt 

Combining (pi5|) with ( pi4|) we get 

J 1 (t)^ 1 - Sign(t) t=( * t<0 . (B16) 
w 2 1 t > v ; 

Thus finally 



= (L ft cos (kL O(fc-i) • (B17) 

V27T2 V 4/ 

This recovers the semiclassical contribution of the bouncing balls 0. There are no 
diffraction corrections to d bhl (k) that are larger than the semiclassical error of the 
bouncing ball contributions. We emphasize that the asymptotic approximation of F(t) 
was invoked after the integration was performed. This is completely justified, because 
the argument of both limits is 0(kz) ^> 1. Replacing the integrand with the asymptotic 
form would be unjustified, since the interesting region near the circle is not asymptotic. 
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For the double repetition the situation is more interesting. To leading order in k we 



have 

^bb,2(^) = 

The relevant integral is 



3 

7T2 



9 



e 2ifcL+if 




dx 



h{t) = 



F 2 (t')dt' = tF 2 (t) 
- sign(t) 



2i ~ . . ! 7T J.2 1 

-F^e 1 ^ - 
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-F(V2t) 



V2tt 

t + V~i/{2n) t<0 
-Vi/(27r) t > 



(B18) 

(B19) 
(B20) 



where the last line was obtained using the asymptotic approximation of F(t). Equation 
(B20) indicates, that in the lit region (t < 0, \t\ ^> 1) we obtain a linear contribution, 
together with a constant that comes from the integration near t ~ 0. In the shadowed 
region we have only a constant contribution from the t ~ region. Inserting (B20) into 
( [BT8D we get 



^bb,2(^) 



3 

7T2 



R cos 2kL 



L 

2r? 



cos(2kL) + Q(k 



(B21) 



Thus, in addition to the semiclassical contribution due to the bouncing balls, we obtain a 
genuine diffractive contribution. It is 0(k°) which is the same as for an unstable periodic 
orbit and therefore must be retained. The explicit form of hit) in (B20) indicates that 
indeed this diffraction contribution is obtained from the region near the circle and is 
not a result of a more global effect. 
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